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abstract 

We present an application of the thermal Comptonization/disk reprocessing model recently proposed by 
Zdziarski, Lubinski, & Smith. We show that the absence of strong optical variations in the presence of strong 
concurrent X-ray variations, similar to those found by HST/RXTE monitoring observations of NGC 3516, can 
be explained by changing the geometry of the Comptonizing plasma rather than the accretion disk itself. The 
total X-ray luminosity of the Comptonizing plasma must decrease as its spatial extent increases. In contrast, the 
disk inner radius must be roughly fixed in order not to produce optical/ultraviolet color variations stronger than 
observed. By including emission due to internal viscous dissipation in the disk, we can roughly match the optical 
and X-ray flux levels and variability amplitudes seen from NGC 3516 during the HST/RXTE campaign. 

Subject headings: galaxies: active — galaxies: individual (NGC 3516) — galaxies: Seyfert — X-rays: galaxies 



1. INTRODUCTION 

A ubiquitous feature of the optical/UV emission from Type 1 
Seyfert galaxies (Syls) is the highly correlated variability of the 
continuum flux across a wide range of wavelengths (e.g., Peter- 
son et al. 1998). When this behavior was first observed during 
the 1989-1990 broad line region reverberation mapping cam- 
paign for NGC 5548 (Peterson et al. 1991; Clavel et al. 1991), 
it was immediately recognized that the upper limits for the rela- 
tive lags between the optical and UV wavebands were too small 
for the variability to have been mediated by hydrodynamic pro- 
cesses in the putative accretion disk — the differences in the 
characteristic radii for the emission in the respective wavebands 
were too large. In fact, the upper limits for the interband lags 
of < 2 days required signal speeds of > 0.1c which could not 
have occurred within the disk itself. This fact and the obser- 
vation that the optical/UV continuum became bluer as the flux 
increased led to the conclusion that the continuum variations 
were due to thermal reprocessing in the disk of higher energy 
radiation produced by a single varying source, such as the X-ray 
source posited to exist at the disk's inner regions, near the cen- 
tral black hole (KroUk et al. 1991; Courvoisier & Clavel 1991; 
Collin-Souffrin 1991). Subsequent observations of other Syls 
showed that highly correlated broad band optical/UV variabil- 
ity is generic to these objects (e.g., Courvoisier & Clavel 1991). 

Recently, the thermal reprocessing picture has been chal- 
lenged by simultaneous optical, UV, and X-ray observations 
of Syls. The first object so observed was NGC 7469 by the 
lUE and RXTE satelHtes over ~ 30 days in 1997 (Nandra et 
al. 1998). Although the UV and X-ray light curves showed 
some similar variability characteristics both in amplitude and 
time scales, these light curves were not nearly as well correlated 
as would be expected if the UV emission was due to thermal 
reprocessing. Naively, the variations in the UV should follow 
those in the X-rays with some characteristic time delay, but the 
maxima in the UV light curves preceded similarly shaped max- 
ima in the X-ray light curve by ^ 4 days, while minima in both 
light curves occurred nearly simultaneously. 

The 1998 observations of NGC 3516 by HST, ASCA, and 



RXTE seem to be equally troubling for the thermal reprocess- 
ing model (Edelson et al. 2000). These observations consisted 
of intensive monitoring in optical, UV, and X-ray bands over 
3 days and so were able to probe correlations on far shorter 
time scales than any previous set of observations. The X-rays 
showed very strong variations (^ 60% peak-to-peak), while the 
changes in the optical continuum were much smaller (^ 3%). 
Because of the accurate relative photometry afforded by the 
HST/STIS CCDs, the measurement uncertainties were suffi- 
ciently small to reveal significant variability that was strongly 
correlated across the optical bands. Hence, the original motiva- 
tion for the thermal reprocessing model was still present, but as 
for NGC 7469, the X-ray and optical flux variations were not 
strongly correlated on any relevant time scale. 

More recently, Nandra et al. (2000) have re-examined the 
1997 observations of NGC 7469 and found that while the rela- 
tively narrow band 2-lOkeV X-ray flux wasn't well-correlated 
with the UV flux the spectral index of the X-ray emission was. 
This led them to suggest that the UV emission is actually corre- 
lated with the bolometric X-ray flux. In this Letter, we expand 
upon that idea and describe what is required for the multiwave- 
band phenomenology exhibited by NGC 35 16 to be due to disk 
thermal reprocessing in the context of a thermal Comptoniza- 
tion model for the X-ray emission. 

2. X-RAY THERMAL REPROCESSING IN THE SPHEREH-DISK 
GEOMETRY 

Most of the important details of the model are described in 
Zdziarski, Lubinski, & Smith (1999; hereafter ZLS) which the 
interested reader is encouraged to consult. We will proceed by 
summarizing the significant features and results from that work. 
This model was proposed by ZLS to explain the apparent cor- 
relation between the X-ray spectral index, F, and the relative 
magnitude, R, of the so-called Compton reflection hump seen 
in the hard X-ray spectra of Syls and X-ray binaries. ZLS 
argued that this correlation results from the reprocessing and 
Comptonizing interactions between the disk and the hot plasma 
believed to produce the X-ray continuum. Expanding upon the 
work of previous authors (e.g., Poutanen, Krolik, & Ryde 1997; 
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Zdziarski et al. 1998), ZLS noted that the R-T correlation could 
arise in a geometry in which the Comptonizing plasma occu- 
pies a spherical region of radius r^ centered on the black hole. 
The accretion disk has an inner radius r^in which penetrates this 
sphere at its equator by varying amounts. This varying penetra- 
tion could be due to changes in r,, r^in, or both. Smaller val- 
ues of r^in/rs result in a larger fraction of the X-ray photons 
from the sphere being intercepted by the disk. The covering 
factor, R, for Compton reflection and thermal reprocessing is 
larger, and a larger fraction of the disk soft photons re-enter the 
plasma precipitating additional cooling via Compton losses and 
thereby reducing the temperature of the plasma. This results in 
a softer X-ray spectrum and a larger value for the photon index 
r. Enforcing energy balance under the assumption that the disk 
emits only reprocessed radiation, ZLS derived a R-T relation 
which looks very similar to the observed correlations. 

We follow ZLS by assuming that the region of Comptoniz- 
ing plasma is homogeneous and optically thin and that we can 
neglect internal viscous dissipation in the disk, at least initially. 
Reprocessing of incident X-rays by a razor-thin disk then gives 
a radial effective temperature distribution 

-,1/4 



magnitude of the X-ray spectrum, we use another empirical re- 
sult of Beloborodov (1999a) to estimate the temperature of the 
Comptonizing plasma, F « (9/4)y"^''^, where the Compton pa- 
rameter is given hy y = 4(0(. + 40^)t(t+ 1) and 9e = ksTe/meC^. 
We assume a typical value for the Thomson depth, t = 1, and 
solve the preceding equations to find the electron temperature 
Te. We model the X-ray continuum as a power-law with an ex- 
ponential roll-over, Lg oc (E/Ecy~^ exp(-E /Ec), where the cut- 
off energy is taken to be Ec = 2kBTe appropriate for a Thomson 
depth T ~ 1 . The normalization, Lq, of the X-ray spectrum is 
given by 



Lv = Ln 



(E/E,y-' exp(-E/E,)d(E/E,). 



(3) 
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For F < 2, the precise value of the minimum photon energy 
is not very important so we have set it at a nominal value of 
£ . =0.01keV. 



T(r) = 



-h 



(1-fl) 



(1) 



Here L^ is the total luminosity of the Comptonizing plasma, ag 
is the Stefan-Boltzmann constant, and a is the disk albedo. The 
dimensionless function h{r/rs) is given by Eq. 1 of ZLS. 

Assuming local blackbody emission and that the optical por- 
tion of the observed spectrum originates at large radii, we may 
use the asymptotic form h{r/rs S> 1) ^ {TT/A){rs/ff' in Eq. 1 to 
obtain the disk spectrum 
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3. COMPARISONS TO OBSERVATIONS 

We have computed thermal Comptonization and disk repro- 
cessing spectra in this model for NGC 3516. The X-ray lu- 
minosity L^ is found by fitting the disk spectrum to the mean 
optical fluxes measured by HST (Edelson et al. 2000). We fix 
the disk inner radius at rmin = 3 x lO'^ cm and plot the result- 
ing spectral energy distributions (SEDs) in Figure 1 for values 
of the sphere radius in the range r, = 1.5-5 x lO'^cm (solid 
curves). As claimed, the calculated disk spectrum is strikingly 
constant, with only relatively small variations in the bluer parts 
of the spectrum. For comparison, we plot the mean fluxes mea- 
sured by HST. Given the crudeness of this model, the quality of 
the fit is satisfactory even for the UV flux. ' 

Since we have thus far neglected local viscous dissipation in 
the disk, the black hole mass does not enter into the calcula- 
tion of the disk emission. However, the value of the disk inner 
radius can, in principle, be constrained by the optical and UV 
data. The dashed curve in Figure 1 is the SED obtained for 
''min = Tv = 6 X lO'^ cm. This shows that while the disk spec- 
trum is relatively insensitive to factor ^ 3 changes in r, for a 
given optical flux, comparably sized changes in r^in produce 
strong color variations in the UV. Therefore, in order to ac- 
count for the UV flux, the disk inner radius cannot be too much 
larger than rmin = 3 x 10'^ cm. For a central black hole mass 
of 3 X 10^ M0(e.g., Edelson & Nandra 1999), fliis value for the 
disk inner radius is roughly consistent with innermost radius 
implied by ASCA observations of the broad Fe Ka emission 
line (Nandra et al. 1999). 

Despite the fact that the optical/UV continuum barely 
changes as a function of r,, the X-ray spectrum varies dras- 
tically. In the 2-lOkeV band in particular, the flux changes 
by factors of ~ 2 for the spectra shown. To relate this result 
more concretely to the X-ray observations of NGC 3516, we 
plot in Figure 2 the X-ray spectral index versus 2-lOkeV flux 
for a family of curves corresponding to different optical fluxes 
with the position on any given curve being parameterized by the 
value of r,. We propose that the HST/ASCA/RXTE monitoring 
observations of NGC 3516 can be explained in this model by 
changes in r, and accompanying changes in L^ so that L^ oc r~^ 
is approximately maintained, and the source takes a path in the 
r-^2-10 plane which lies nearly along a curve of constant opti- 
cal flux. Any significant variations in the optical/UV continuum 

' This datum was not included in the fit for Lx since the UV flux was only monitored over a significantly shorter period than the optical data due to SAA passages. 
In addition, the UV data point, at ~ 1350A, is in a part of the spectrum which may be contaminated by the broad wings of the SilV emission line (see Edelson et al. 
2000, Fig. 3). Unfortunately, similar contamination problems may exist for the optical data as well (see Edelson et al. 2000, Fig. 4). 



exp{E^ryykBA) 

where T = Ar^'^, A = [3Lj,r,{l-a)/64TTaB]^/*, and E^ is the 
photon energy. Neglecting variations in the disk albedo and as- 
suming that rmin is fixed on the time scales of interest, a constant 
optical flux implies a constant value of A. In the presence of 
variations of the X-ray luminosity, this requires Ljc oc rj'. 

For r 2± Tj, i.e., at higher frequencies, it turns out that the large 
radius limit of h{r/rs) is a surprisingly good approximation for 
evaluating the disk spectrum since it enters as the 1/4-power in 
r, and at r/r^ = 1, the function takes the value h{\) = 4/3 (ZLS), 
and we have (7r/4)i/'* « (4/3)'/'* « 1. Hence, for r, ~ r„i„, 
adopting the asymptotic form for liir/r,) throughout and set- 
ting Li ex r~^, all rj-dependence drops out of Eq. 2, and in this 
approximation, the resulting disk spectrum has a constant shape 
and magnitude for different values of r,. The exact expression 
for h{r/rs) flattens for r < rs, so when rmin < fs, the disk spec- 
trum predicted by the full calculation will be somewhat softer, 
but not appreciably at the wavelengths which have actually been 
observed. 

Given the X-ray luminosity L^, we estimate the X-ray spec- 
trum produced by thermal Comptonization in the sphere. For 
this we use a second result from ZLS which relates the am- 
plification factor for Comptonization, A, to the ratio rmin/r? 
(ZLS, Eq. 2). Following ZLS, we use the empirical result of 
Beloborodov (1999a) which gives the X-ray spectral index as 
a function of A: F « 2.33(1 -A)"'-"". To determine the overall 



will correspond to deviations from the constant flux curve and 
therefore from the L^ oc rj^ relation. Furthermore, the changes 
associated with such departures will be broad band and will 
show the same color variations that have been previously asso- 
ciated with thermal reprocessing. 

For comparison, we have plotted in Figure 2 values of F and 
F2-10 which we have obtained from fitting the publicly archived 
RXTE NGC 3516 data for the observations in question. Al- 
though this simple model overpredicts the X-ray flux by a factor 
^ 4 for the observed range of spectral indices, the correlation 
of spectral index and flux given by the model is strikingly sim- 
ilar to that of the data. The flux discrepancy could be reduced 
somewhat by selecting a smaller value for the Thomson depth 
T since for a given spectral index (and y-parameter) that will 
imply a larger value of Tg and more of the X-ray luminosity can 
be "hidden" at higher energies. Alternatively, a large reduction 
can be obtained by including contributions to the disk emission 
from local viscous dissipation. This corresponds to adding a 
term ex (GMM / r^)(l -(ri/ry^^) inside the square brackets in 
the expression for the disk temperature (Eq. 1) where r/ = r^ for 
a maximal Kerr metric and rj = 6rg for a Schwarzschild met- 
ric. Since the viscous dissipation term is independent of r,, it 
acts essentially as a constant offset for the required value of the 
X-ray luminosity, and L^ ex rj' is still needed to maintain a con- 
stant optical flux. In order to match the observed X-ray fluxes, 
roughly 1/4 of the total disk luminosity must be due to thermal 
reprocessing with the remaining 3/4 being due to viscous dissi- 
pation. Assuming M = 3 x 10^ Mq and r/ = rg, we have solved 
for the energy balance condition and find that an accretion rate 
of M ~ 3.6 X 10"^ M0 yr"' is required to match the mean X-ray 
flux. This solution is plotted as the thick solid curve in Figure 2. 

Although the energy balance condition we have used only 
restricts the properties of the Comptonizing plasma, global en- 
ergy balance is roughly satisfied by this solution. The bolo- 
metric X-ray luminosity which gives the observed mean X- 
ray flux and spectral index is L^ ~ 7.0 x 10"*-' erg s"', while 
the above accretion rate implies a disk luminosity of Ldisk — 
3.3 X 10"*^ erg s"' due solely to viscous dissipation. For a disk 
extending all the way to rj = rg (rather than being truncated at 
r,), the total implied luminosity is Ltot — 1.0 x lO'^'^erg s"', and 
we have 

hot ^Lx+ Ldisk- (4) 

This would be expected if the accretion luminosity that would 
otherwise emerge as local black body emission from the inner 
regions of the untruncated disk were instead partitioned to the 
heating of the Comptonizing plasma. We note that we have 
arbitrarily set the Thomson depth to r = 1 in computing this so- 
lution. By fine tuning the values of t and M, we find that we 
can satisfy the global energy balance condition arbitrarily well 
while simultaneously matching the observed mean X-ray flux, 
spectral index and optical continuum levels. 

4. DISCUSSION AND CONCLUSIONS 

The fact that an approximately constant optical/UV contin- 
uum can in principle result from thermal reprocessing of a 
highly variable X-ray flux is connected with the fact that the 
radial distribution of the flux incident upon the disk is well ap- 
proximated by a power-law, even very close to the innermost 
disk radius. As we showed above, the constant flux constraint, 
in conjunction with the assumed constancy of the disk inner 
radius and disk albedo, causes the radius of the Comptoniz- 
ing sphere to drop out of the integral describing the disk spec- 



trum. However, this constraint does not demand that the shape 
of the incident flux distribution be a power-law. The only re- 
quirement is that the relevant physical parameters affecting the 
incident flux can be considered independently from the expres- 
sion which describes the shape of the distribution. For example, 
consider the model proposed by Beloborodov (1999b) in which 
the Comptonizing medium consists of a blob of plasma which 
is at some height r/, above the disk. The blob moves with mildly 
relativistic velocities either vertically away from or towards the 
underlying accretion disk, and the effective reflection fraction 
is controlled by relativistic aberration of the flux from the blob. 
Using the formulae in Beloborodov (1999b), the radial distri- 
bution of the flux incident upon the disk is given by 






(r^+1) 
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/3 + (r2+l)2 
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where /? is the blob velocity in units of c, 7 = (1 - Z?^)"'/^, and 
f = r/rh is the disk radius in units of the blob height. For the 
range of blob velocities required by this model to reproduce 
the observed R-T correlation, /3 ~ 0.3-0.7, the factor in the 
square brackets which gives the shape of the flux distribution 
is largely insensitive to the exact value of (3. Therefore, since 
the lower limit of the radial flux integral can be taken to be zero 
(cf. Eq. 2), a nearly constant optical flux in this model can be 
obtained if L^ « r\"i'^. 

Although neither of the aforementioned models offer any 
reason for the luminosity scalings which seem to be required by 
the NGC 3516 observations, any plausible explanation would 
still have to contend with why NGC 7469 and probably other 
Syls violate it. Therefore, we will not attempt to posit a phys- 
ically motivated rationale for the specific relationship between 
X-ray luminosity and geometry required for any particular ob- 
ject, but a cursory examination of the important time scales 
may provide some insights. For NGC 3516, the typical vari- 
ability time scale exhibited by the 2-10 keV light curve is 
^ 20-30 ks (Edelson et al. 2000). Assuming a radial Thom- 
son depth r = 1 and r, ~ rmin, the cooling time is very short: 
fcooi ^ NekBTe/L, < 10 s, where A^^ = (4TrPj3)(T/r,aT) is flie 
total number of electrons in the sphere with mean tempera- 
ture Te ^ lO^K. Unfortunately, the thermal expansion time for 
a sphere of this size is fairly long, approximately fexp ^ 100 ks. 
This assumes that the protons and the electrons are at the same 
temperature. If we have a two-temperature plasma in which the 
protons possess a virial temperature of ^ 10'^ K, then the ex- 
pansion time scale is reduced substantially. Alternatively, this 
could also be achieved if we have a pair plasma instead. 

We have assumed for simplicity that the thermal emission 
from the disk is locally blackbody. Detailed modeling of bare 
disk atmospheres appropriate for Syls indicates that this is a 
moderately reasonable approximation in the optical, but signif- 
icant deviations can occur at higher frequencies (e.g., Hubeny 
et al. 2000, 2001). This appears to remain true for disks pow- 
ered exclusively by external X-ray illumination (Sincell & Kro- 
lik 1997). Our treatment of the reprocessed flux in the optical 
is therefore probably fairly good, but the full seed photon spec- 
trum illuminating the Comptonizing plasma will not be very ac- 
curate. More detailed modeling of the disk spectrum will likely 
affect the specific values we deduce for the disk inner radius. 
X-ray luminosity, and Comptonizing sphere radius, but it will 
not change our general conclusions. 

The existing simultaneous optical, UV, and X-ray data plus 
the global energy balance condition (Eq. 4) provide just enough 



information to infer unique values for the five model param- 
eters, Tmin, rs, Te, T, and M, required to describe the entire 
spectral energy distribution. Therefore, a crucial test of this 
general scheme for linking thermal Comptonization and ther- 
mal reprocessing will be provided by simultaneous optical, 
UV, and broad band X-ray monitoring observations which in- 
clude reliable measurements of the thermal roll-over at ener- 
gies > 100 keV. These latter measurements are necessary for 
making estimates of the X-ray luminosity L^ which can then be 
compared with the luminosities required by the model. Forth- 
coming INTEGRAL observations of Syls will provide some of 
the necessary spectral coverage to address this issue, but what is 
really required is a hard X-ray/soft gamma-ray telescope which 



can measure the X-ray spectrum and thermal roll-over on time 
scales comparable to the shortest variability time scales seen 
from these objects. 
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Fig. 1. — Model SEDs which have been fit to the HST optical fluxes for NGC 3516. For the solid curves, the disk inner radius 
is set at r^[^ = 3 X lO" cm, the Thomson depth of the Comptonizing sphere is t = 1, and the radius of the sphere takes values 
r, = 1.5-5 X lO'^ cm. The smaller values of r, correspond to the harder X-ray spectra. The dashed curve is the SED corresponding to 
''mill = Tj = 6 X lO" cm. The dotted lines indicate the 2-10 keV band. 
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Fig. 2. — X-ray spectral index versus 2-lOkeV flux for families of curves each of which produces a different constant optical 
continuum level. The rightmost curve represents the model calculations shown in Figure 1 and the next four curves to its left have been 
computed for 87.5, 75, 62.5, and 50% of the mean optical flux level for NGC 3516. The X-ray emission can vary significantly and still 
produce a nearly constant opticalrtJV continuum if the source follows a path along one of these curves such as the one labeled "NGC 
3516". The thick solid curve has been fit to the full optical flux, but emission due to internal viscous dissipation in the disk has been 
included to reduce the required X-ray luminosity. The neighboring dashed curves bound the region of ±1.5% optical variations. The 
data points with error bars are from spectral fits to the archived 1998 RXTE data for NGC 3516. 



